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The derivative discontinuity in the exact exchange-correlation potential of ensemble Density Func- 
tional Theory (DFT) is investigated at the specific integer number that corresponds to the maximum 
number of bound electrons, J m ax- A recently developed complex-scaled analog of DFT is extended 
to fractional particle numbers and used to study ensembles of both bound and metastable states. 
It is found that the exact exchange-correlation potential experiences discontinuous jumps at integer 
particle numbers including J max - For integers below Jmax the jump is purely real because of the real 
shift in the chemical potential. At J m ax, the jump has a non-zero imaginary component reflecting 
the finite lifetime of the {Jmax + 1) state. 



INTRODUCTION 

The ground-state energy of an open system as a func- 
tion of particle number has derivative discontinuities at 
integer values and is linear in between integers [IHSj (see 
Fig. [I). Since Density Functional Theory (DFT) [1(5] is 
formally exact, an extension of the theory to fractional 
particle numbers should also include this behavior. Most 
approximate functionals used in practice have non-linear 
behavior between integers and lack derivative discontinu- 
ities in the energy. These incorrect features of approxi- 
mate functionals have profound consequences [3] • The as- 
sociated derealization and static correlation errors typ- 
ically lead to an underestimation of chemical barriers, 
band gaps, long-range charge transfer excitations, and 
energies of dissociating molecules and an overestimation 
of binding energies. 

The derivative discontinuities in the energy have been 
studied at integer particle numbers lower than the max- 
imum number of bound electrons. For example, by 
considering the 2-component ensemble formed by a J- 
electron system and the corresponding (J -I- l)-electron 
system in an exactly-solvable 2-electron model system, 
Sagvolden and Perdew established that the exact Kohn- 
Sham exchange-correlation potential experiences a con- 
stant positive jump at integer J = 1 !7|. They also pos- 
tulated that one could study the discontinuity at the in- 
teger number, J max , of particles that corresponds to the 
maximum number of bound electrons. However, their 
useful analysis could not be extended to J max because 
the {Jmax + 1) electron system is not bound. The discon- 
tinuity, on the other hand, is a property of the potential 
at the integer value where the system is bound, so it 
would be useful to study such behavior. 
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FIG. 1. Cartoon illustrating the ground-state energy of an 
open system as a function of particle number. If a long-lived 
metastable state exists above the maximum number of elec- 
trons {Jmax), the curve can be thought to turn upward be- 
cause of the negative electron affinity of the Jma^-electron 
system. 



The electron affinity of a J-electron system is defined 
as the energy difference, A = Ej — Ej + i, when one elec- 
tron is added to the system. In the case when J = J max , 
the affinity could be taken as zero because the {J m ax + 1) 
system is not stable (dashed line in Figjl]). However, the 
{Jmax + 1) system could form a long-lived metastable 
anion, and such a system is often described as having a 
"negative electron affinity" (NEA) (HE]. Therefore if the 
energy is plotted as a function of particle number (TV), 
below Jmax we have a concave up curve with discontinu- 
ities at the integer values ( J's) and negative slopes, but 
beyond J m ax the curve can be thought to turn upward 
because of the change in sign of the electron affinity (see 
Fig [I]). One physical system that exhibits this behavior 
is the nitrogen molecule {J m ax = 14) whose negative ion 
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is unbound but long-lived. The NEA of N 2 is -2.3 eV 
or -0.084 Hartree because of the existence of a long-lived 
(14+l)-electron metastable state (N%) 2.3 eV above the 
ground state of N 2 . The width of this N% state is 0.57 
eV corresponding to a lifetime of 1.2 femtoseconds. 

When a molecule with J max electrons is in contact with 
an infinite but distant metallic reservoir of work func- 
tion W smaller than the magnitude of the NEA of the 
molecule (W < \A\, e.g. N 2 in contact with a slab of 
Cs), the ground-state density of the combined system is 
particularly interesting. At infinite separation, it is sim- 
ply the sum of neutral, isolated N 2 and metallic densi- 
ties. At finite separation, however, close to the molecule, 
it is given by a linear combination of the J ma:E -electron 
ground-state density (for N 2 ) and the ( J max + l)-electron 
density of metastable . 

As mentioned, the lack of the correct derivative dis- 
continuities in approximate functionals has a dramatic 
impact on the calculation of electronic properties, and it 
is unknown how the presence of a long-lived metastable 
state influences the exact exchange-correlation potential 
as the particle number goes through J m ax- Although the 
accurate calculation of electron affinities has been ex- 
plored with approximate DFT [TO], it has been unclear 
if one should use negative electron affinities (as opposed 
to setting the affinity to zero) in the calculation of prop- 
erties such as the chemical hardness which depend on 
derivatives of the energy with respect to particle num- 
ber [TT]. If the behavior of the energy functional was 
known at fractional particle numbers above J ma x, one 
could formulate answers to such questions. 

In this work, we use complex "densities" calcu- 
lated from complex-scaling theory [T2"HT4"] to study the 
derivative discontinuity in the Kohn-Sham exchange- 
correlation potential at the maximum number of bound 
electrons. Since the complex energy functionals used give 
both the bound energies and metastable {J ma x + 1) en- 
ergies and lifetimes, the method presented here can be 
used to probe the behavior of the groud-state exchange- 
correlation potential around J max . 

In previous work, we have demonstrated that the en- 
ergy and lifetime of the lowest metastable state can be 
extracted from a complex density with a properly scaled 
energy functional [THIHH]. Also, a complex "Kohn-Sham" 
system can be constructed that facilitates self-consistent 
calculations on many-electron systems with the complex 
density as the primary variable |17| . Zhang et al. have 
developed a Levy-Nagy extension of this formalism to 
treat higher energy resonances (or excited states of the 
metastable system) |18j . and in related work Maitra et 
al. have considered autoionizing resonances within time- 
dependent DFT [19]. Also, Ernzerhof and co-workers 
have developed an approach applicable to molecules con- 
nected to metallic leads where complex absorbing poten- 
tials are added within a complex-DFT framework [2"Ull21| . 
However, the complex potentials in the "Density Func- 



tional Resonance Theory" (DFRT) of Ref. [T7] result 
from a variational calculation, and they are obtained self- 
consistently for the TV-electron system treated as isolated, 
rather than added to the Hamiltonian from the start to 
model an open system. 



DENSITY FUNCTIONAL RESONANCE 
THEORY AND COMPLEX DENSITIES 

We use the complex analog of Kohn-Sham DFT, Den- 
sity Functional Resonance Theory (DFRT), to treat 
metastable systems [TT]. The lowest-energy resonance 
energy and lifetime of a system is encoded into a corre- 
sponding complex resonance density [15j defined by, 

n e (r) = <*£|n(r)|*£) (1) 

where n(r) is the density operator, and (^Sg\ and 1^^) 
are the left and right eigenvectors of the complex-scaled 
Hamiltonian, Hg, corresponding to the Lowest-Energy 
Resonance, or LER (see Ref. [22] for a review of complex- 
scaling theory) . The angle 9 is the complex-scaling angle 
in the transformation of the coordinates from r to re . 
We require that ng(r) be normalized to the number of 
electrons, as real densities are: 

J dr n 8 (r) = J (2) 

The lifetime, £, of the resonance is defined as 
(2 Im(Ee))^ 1 , where Eg is a complex eigenvalue of Hg 
corresponding to a pole in the scattering matrix. The 
real part of Eg, £, will be referred to as the "resonance 
energy." 

There is a one-to-one correspondance between 
complex-scaled external potentials and their correspond- 
ing LER complex densities [16J. Therefore, the energy 
and lifetime of the LER can be extracted from ng with 
a properly scaled energy functional. For J electrons this 
functional is 

£[ng] - % -L-\ng\ = T e s [n e ] + J dr n g (r)v(re ie ) 

+E°[n 9 ]+E° [ne] (3) 

in analogy to standard KS-DFT, with T®[n e ] = 
e- 2ie T s [ng] and E^[n e ] = e- iB E n [n e ], where T s [n e ] and 
E H [ng] are the standard non-interacting kinetic energy 
and Hartree functionals evaluated at the complex densi- 
ties. Eq. [3] then defines E^ c [ng}. 

The system of interacting electrons whose LER den- 
sity is ng(r) is mapped to one of J particles mov- 
ing independently in a complex "Kohn-Sham" potential 
v s (r) defined such that its J occupied complex orbitals 
{4>®(r)} yield the interacting LER-density via ng(r) — 
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sent at ion 



L \n(r)\^ R 



In Moiseyev's Hermitian repre- 
, the complex Kohn-Sham equations are: 



hi - e t -h 2 - 2TT 1 \ / Re(^) 
12 + 2rr l k - Ei )\ Im(^) 



(4) 



where hi = - \ cos(26»)V 2 + Re(uf(r)), and /i 2 = 
| sin(26»)V 2 + Im(uf (r)). The set of {ej and {t,} pro- 
vide the orbital resonance energies and lifetimes of the 
Kohn-Sham particles. 

The complex variational principle |22j along with the 
assumption that the orbitals used to construct the den- 
sity can be expanded in an orthonormal basis leads to 
the Euler-Lagrange equation: 



SE g [n e ] 
Sn g 



H J drng(r) = 



(5) 



Performing the variation in Eq. [3] and comparing with 
Eq.|4]leads to an expression for the Kohn-Sham potential 
that is again analogous to that of standard KS-DFT: 

v e s (r) = v(re w ) + e~ ie v n [n e ](r) + v e xc [n e ](r) , (6) 

where v^ c [n e ](r) = SE^ c [ne}/5n e (r)\ LER . 

DFRT FOR FRACTIONAL PARTICLE NUMBERS 

We start by considering the Levy-Lieb constained- 
search functional that has already been applied to both 
pure and ensemble states. 



min (T + V, 



ee/r 



(7) 



where T is the kinetic energy operator, V ee is the electron- 
electron repulsion, and T is a statistical mixture or en- 
semble of pure states, T — J2i \^i)Pi(^i\- The sum of 
the probabilities Pi is defined to be 1, and the expecta- 
tion value of some observable is: 



(6>=^p i <* i |6|* i >=Tr 



TO 
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For a fractional particle number N between J and ( J+ 
1), the ensemble of a J-electron pure state and a ( J + 1)- 
electron pure state minimizes Eq. [7] This is true as long 
as the energy as a fraction of particle number obeys the 
"concave-upward" condition, or Ej < (Ej_i + E J+ i)/2, 
for which no counter-example has been found in nature 
(but the condition has not been proven in general). 

In extending DFRT to fractional particle numbers, we 
define an ensemble of pure complex- scaled states as: 



(9) 



where has the same meaning and follows the same 
conditions as the usual, non-complex-scaled, case. Then, 
the bi-expectation value is defined as 



On 



T e O 



(10) 



We now make the assumption that the resonance of 
lowest energy is also the one with the longest lifetime. 
Like the convexity assumption above, this has never been 
proven, but it has been found empirically to be the typi- 
cal case [24. 25j. The constained-search functional is then 
defined by: 



F^ h [n$] = min 



Re 
-21m 



Tr 



T e (e~^T + V: e ) (11) 



where V® e is the complex-scaled electron-electron inter- 
action operator. The ground-ensemble energy is: 



E e [n e ] = 



mm 



Re 
-21m 



F?M + / dv n e (r)v(re id ) 



(12) 

If the specific statistical mixture is that of the J-clcctron 
state, *&gj, and the (J + l)-electron state, ^e.j+ii the 
complex ground-ensemble energy for a fractional particle 
number N ( J < N < J + 1) is, 

E e (N) = [1 — (JV — J)]E e j + (N - J)E e }+1 (13) 

where Ej is the bound ground state (or if no bound state 
exists, lowest metastable state) complex energy of the J- 
electron system, and likewise Ej , 1 is the bound ground 
state (or if no bound state exists, lowest or second-lowest 
metastable) complex energy of the J + 1 electron sys- 
tem. Again, the lifetime, £, of the state is given by 
(2 Im(Eg))^ 1 (this will be zero if the state is bound), 
and the real part of Eg, £, is the ground-state energy if 
the system is bound or the "resonance energy" if the sys- 
tem is unbound. The complex ground-ensemble density 



<(r) = [1— (N — J)} n J e (r) + (N — J)n^(r) (14) 
One might ask whether the ensemble of the J and (J + l) 



electron states minimizes Eq. 12 or whether we have a 
"concave-upward" condition given that we now are deal- 
ing with complex energies. The real part of the energy 
(position of the resonance or bound state) follows the 
concave-upward condition already discussed. However, 
assuming that the LER has the longest lifetime, the imag- 
inary part of the energy as a function of particle number 
would be concave-down (at least up to J max + 1). Yet, 
the constrained search in Eq. [12] is defined as giving the 
minimum width, T, or the maximum lifetime, and the 
width docs follow a "concave-upward" condition. 

This extension of DFRT allows us to study the deriva- 
tive discontinuity at integer values below J m ax, because 
DFRT can be applied to bound states, but it also al- 
lows us to investigate the discontinuity (if there is one) 

&t Jmax- 
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SIGNIFICANCE OF THE COMPLEX 
KOHN-SHAM ORBITAL ENERGIES 



system with one decay channel discussed above, compar- 
ing with Eq. |16[ 



One related topic that needs to be discussed before 
studying the discontinuity at J ma x is the physical mean- 
ing of the Kohn-Sham orbital energies. One can do this 
by looking at the assymptotic behavior of densities. If 
we consider the tails of a real, bound ground-state den- 
sity, we know that they decay exponentially like e~ 2 ^ r , 
where I is the ionization potential [27] . The highest 
occupied Kohn-Sham orbital energy is then given by the 
negative of /. However, in DFRT we must take care- 
ful note of how a resonance wavefunction, cf>, decays. As 
r — > oo, 



Ce 1 



where, 



a ,3 



2(E a - Ef) 



And the complex-scaled density would decay as: 



(,;) 



C 2 



(15) 



(16) 



(17) 



In these expressions the index a defines the eigenvalue 
of the complex-rotated Hamiltonian, j labels the decay 
channel, and C is a constant [25]. Ej h is the thresh- 
old energy of the decay channel, and E a is the complex 
eigenvalue of the resonance. The wavefunction decay is 
governed by the energy difference (E a — E* h ). 

For simplicity of discussion, consider the LER of a 
( J + l)-electron unbound system that decays to a bound 
J-electron system (this is the single dominant channel 
of decay) . £f would equal the energy of the J-electron 
state (Ej). Therefore the tail of the complex density 
of the resonance behaves like exp (i2W2(Erj + u — Ej)r) , 
where Erj+i) = £u+i) — ?i is the complex LER reso- 
nance eigenvalue {£ (j + i)-position or energy, T-width or 
inverse lifetime). The energy difference in the exponen- 
tial tail of the complex density can then be related to the 
NEA of the J-electron system, A — (Ej — £/j +1 \): 



r . 



(E a - Ef) = (£ {J+l) - -i - Ej) = {-A - -i) , (18) 

For a ( J + l)-electron unbound system with multiple 
decay channels (or a non-sequential decay process) it is 
useful to look at the time-dependence of the states (par- 
tial widths and branching ratios) to gain understanding 
of the decay process [28]. The relation to a NEA is not 
immediately clear. 

The complex density of the LER in DFRT, constructed 
from the KS orbitals, behaves exactly as in Eq. [17] with 
the wavenumber given in terms of Kohn-Sham quantities: 

k *,3 = \J 2 ( e l - £ f )■ ef is the KS "threshold energy," 

and e a is the complex HOMO resonance energy. For the 



(4 



= (-A-U) 
V 2 ' 



(19) 



This provides some physical interpretation for e^, but we 
need a clear definition of ef . In analogy to ground-state 
DFT, one can write the LER's complex energy in terms 
of a sum of KS complex orbital energies [17] : 



£(.7+1) [ n e\ - 



r[n e ] _ e 



J 

E« 

i=l 



EL 



drvi xc (r)n e (v) (20) 



Where the HOMO energy has been pulled out from the 
eigenvalue sum. We can then write the following expres- 
sion for e„, 



£(.7+1) M - ( r M/ 2 )* - £M 



(21) 



where, 



j 

£4 



El 



\n e 



drv e HXC (r)n e (v) (22) 



Comparing Eq. 21 with Eq. 19 



e th = £ {J+1) [n e ] - (T[n e ]/2)t - £[n g } + (A+ -i) 

= Ej - t[ng] (23) 

Therefore, e th is just the threshold energy in the real 
interacting system, Ej, shifted by £. Now with this defi- 
nition of e th we can relate the HOMO energy of DFRT to 
physical quantities via Eq. [19] (for a system with a single 
dominant decay channel). Also, note that the standard 
Koopmans' theorem for DFT is recovered when the sys- 
tem of interest is bound. 

Table [T] summarizes this analysis. 

MODEL PROBLEM 

To study the behavior of the exchange-correlation po- 
tential at J TXWIX 

we consider a system of two interacting 
electrons moving in a one-dimensional potential which 
can be solved exactly using finite difference methods. We 
study a Hamiltonian where the electrons interact via a 
soft-Coulomb potential of strength A: 



H 



E 

i=l 



Id 2 

2dx 



2 + v(Xi) 



A 



y/l + (Xi - X 2 ) 2 



(24) 



with, v(x) = a 



where a, a, 6, c, and d are constants. Note that the 
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System: 



Decay Channel: 



Bound J-electron No decay —I 

Metastable (J + l)-electron J electron bound state {—A — ^i) + e th 

Metastable (J + l)-electron Multiple Dependent on partial widths / branching ratios 



e|, KS-DFRT HOMO energy 
I, Positive ionization potential of a J-elec tron system 
e th , KS "threshold energy" defined in Eq. 



2:! 



A, Negative electron affinity of a J-electron system 
r, (J + l)-electron resonance width (inverse lifetime) 



TABLE I. Summary of the significance of Kohn-Sham DFRT orbital energies with definitions of the relevant quantities. 



exact complex energies are theoretically independent of 
the scaling parameter 9, but numerical approximation 
introduces some ^-dependence. This dependence dissa- 
pears in the infinite basis-set limit, but it can also be 
addressed with a finite basis set or finite grid by calcu- 
lating "^-trajectories" [22] for an optimum 9. This pro- 
cedure was used, but for the above model problem we 
were able to use enough grid points to extinguish most 
of the (9-dependence of the energies. The ensemble den- 
sity, Eq. [l4j for 1- and 2-electron systems can be inverted 
to construct the exact complex KS potential: 



(25) 



where e H is the highest occupied Kohn-Sham orbital. The 
Hartree contribution to the KS potential can be found 
directly from the density, 



«M0»0 = / dx v ee (x'e i6 ,xe ie )n$(x') 



(26) 



The exact exchange-correlation potential can then be 
found by the difference, 



I i8\ i 8 

—v(xe ) + e H 



(27) 



For the bound case and for particle numbers between 
and 2, Sagvolden and Perdew proved that the kinetic 
and Hartree parts of this potential change continuously 
as the particle number crosses an integer [7]. Note that 
these potentials depend on 9, so we only compare po- 
tentials calculated with the same scaling parameter [T7] . 
For certain values of the parameters this potential can 
be made to support 2 bound states. We show that in 
this case v^ c exibits discontinuous jumps for integers 
J < J max- For a different set of parameters, v(x) can be 
made to have a very weakly bound 1-electron state and 





FIG. 2. The two model potentials used in this study. On the 
left, v(x) with a — 0, a = 9, b — 0.5, c = 0, and d = 0. 
This potential has both a I electron and 2 electron bound 
state (when A = I in v ec ). On the right, v(x) with a = 0.75, 
a = 6, b = 0.05, c = 4, and d = 3. This potential has a very 
weakly bound I electron state and a metastable but long-lived 
2 electron state (when A = I in v ee ). 



a metastable but long-lived 2-electron state. For such a 
system, the maximum number of bound electrons, J ma xi 
is one, and the system has a negative electron affinity, 
A = £j max — £j max +i- Calculations on the bound 1- 
electron and unbound 2-electron states then allows 
to be examined for the integer J = J max - 



RESULTS FOR J < J max 

The model potential v(x) with a = 0, a — 9, b = 0.5, 
c = 0, and d = and v ee with A = 1 supports both 
a 1-electron bound state and a 2-electron bound state. 
The external potential is shown in Fig. [2] (left panel). 
The 1-electron bound energy is E\ = —6.38, and the 2- 
electron bound energy is E 2 = —11.84. Both the 1- and 
2-electron complex densities are found exactly with finite 
differences. The ensemble density is formed according 
and the exact 



to Eq. 
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is calculated for different 
values ofN. This model system is similar to the one 
studied in Ref. In that work, the purely real XC 
potential was shown to experience a jump at the integer 
J = 1. We find a similar result when using our complex- 
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FIG. 3. The real and imaginary parts of « xc above and below 
the integer J = 1 when Jmax — 2 (8 — 0.35). The real part of 
the potential jumps by E2 — 2E\, as in ensemble DFT. The 
imaginary part does not experience any discontinuity (solid 
and dashed lines are on top of each other in the bottom pan- 
els). 



FIG. 4. The real and imaginary parts of d xc above and below 
the integer J = 1 when Jmax = 1 (6 = 0.35). It is clear hear 
that the real part jumps by (Re(i?2) — £a), but we must take 
a closer look at the imaginary part to see the discontinuous 
jump of -r/2 = lm(E 2 ). 



of v „ above and below J„ 



1. The structure of these 



density analysis. The change in the complex density for a 
value of N infinitesimally less than or greater than J = 1 
is negligible, yet the chemical potential, or Kohn-Sham 
HOMO energy (see Table [j}, experiences a purely real 
jump of A/i = E2 — 2Ei on either side of J = 1. Also, 
since there is a smooth change in the complex density, 
the change in the Hartree potential is negligible on either 
side of J = 1. Therefore, the complex XC potential must 
compensate for this shift in the chemical potential by a 
positive jump (see Figure [3). There is no jump in the 
imaginary part of the XC potential because there is no 
jump in the imaginary part of y, or the HOMO energy. In 
other words, both the J = 1 and J = 2 states are bound 
and their energies have an imaginary part of exactly zero. 



RESULTS AT J = J max 

The model potential v{x) with a = 0.75, a = 6, 
b = 0.05, c = 4, and d — 3 and v ee with A = 1 supports 
only one bound state with 1-electron and one metastable 
state with 2 electrons. This potential is shown in Fig. [2] 
(right panel). Its steps mimic the centrifugal barriers 
present in 3D potentials which give rise to shape reso- 
nances. The 1-electron bound energy is E\ = —0.86, and 
the 2-electron complex energy of the metastable state 
is Ei = —0.63 — 0.066«. Both the bound 1-electron and 
metastable 2-electron complex densities are found exactly 
with finite differences. The ensemble density is formed 
and v 6 „ is calculated for different val- 



14 



according to Eq. 
ues of TV. Fig.[^snows both the real and imaginary parts 



potentials is similar to the bound potentials in Fig. [3] ex- 
cept for the dramatic scars in the imaginary part from 
the potential steps which influence the lifetime of the res- 
onance. From Fig. [4] one can see that w xc experiences a 
discontinuous jump in its real part, and if we zoom in on 
w xc a j um P m ^s imaginary part also becomes visible (see 
Fig. [4]) . These jumps are to compensate for the complex 
shift in the chemical potential. The HOMO energy just 
below J m ax = 1 is e Hl x = Ex, but the HOMO energy just 
above J max = 1 is e H , 2 = {-{Ex - £ 2 ) - |i) + e th , where 
F is the width of the 2-elcctron resonance (see Table |!|. 
For this system, we can calculate the exact e th accord- 
ing to Eq. [23} Then e H .2 = -0.17 - 0.15i Therefore, 
A/i = {-{Ex - £2) - §0 + e * h ~Ex = 0.69 - 0.15i. Note 
that here the constant jump in the real part of the XC 
potential is positive, as in previous studies, but the jump 
in the imaginary part is negative because of the negative 
shift in the chemical potential from the existence of a 
long-lived metastable state. 

In the previous work by Sagvolden and Perdew [7], 
it was observed that the exact v KC of ground-state DFT 
slightly above an integer was shifted upward from the w xc 
slightly below an integer. This constant upward shift ex- 
tended out to a radius that depended on N, the fractional 
number of particles, and this radius extends further out 
as N approaches the integer from above. We can see the 
same type of asymptotic behavior in the real part of i> xo 
as a function of the fractional particle number, N (see 
Fig. [5]). As N -> J max , the radius out to which Re(w xc ) 
is shitted upward by a constant goes to infinity. Although 
this radius goes to infinity as N gets infinitesimally close 



7 



> 

CD 
DC 




FIG. 5. Asymptotic behavior of v xc as particle number in- 
creases (6 = 0.35). 



to J m ax from above, there is no contradiction with the 
exact condition stating that the Kohn-Sham potential 
should go to zero at infinity. As stated in Ref. [2j [7] , if 
(N—J) is finite there will always be a radius at which the 
constant shift falls off. For larger values of N the shift is 
clearly not a constant, instead shows a rich structure 
resulting from the transient nature of the (J ma x + 1)- 
electron state. 

When one electron is bound and two electrons are un- 
bound the energy as a function of particle number can 
be represented by two piecewise linear graphs. First, a 
graph that is similar to the cartoon (Fig. [T]) in the intro- 
duction (see Fig[6| shows that the real part of the energy 
of an open system experiences discontinuities at integer 
values, including J ma x- Next, a graph of the imaginary 
part of the energy of an open system shows that there are 
discontinuities at integer values of particle number that 
are greater than or equal to J ma x (see Fig [6| . 



CONCLUSION 

As the particle number of a ground-state open sys- 
tem crosses an integer the exact complex Kohn-Sham 
exchange-correlation potential of DFRT extended to 
fractional particle numbers experiences a discontinuous 
jump. At integers below the maximum number of bound 
electrons this jump is purely real and positive compen- 
sating for a purely real shift in the chemical potential. 
At the specific integer that corresponds to the maximum 
number of bound electrons the jump has a non-zero imag- 
inary component. The real part of the jump is positive 
compensating again for the real shift in the chemical po- 
tential, and the imaginary part of the shift is negative 
due to the existence of a long-lived metastable state. 

Since one can recover ground-state Density Functional 
Theory from DFRT by removing the complex transfor- 
mation, we postulate that the exact exchange-correlation 
potential of DFT extended to fractional particle numbers 
should display a discontinuous jump at J m ax, the maxi- 
mum number of electrons the system can bind. Due to 
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FIG. 6. The real part of an open system's energy as a func- 
tion of particle number experiences discontinuities at integer 
values of particle number including J max . The imaginary part 
experiences discontinuities at integer values of particle num- 
ber greater than of equal to J max . This example is for the 
model system with Jmax = 1. (I is the ionization potential 
of the 1 electron system, A is the affinity of the 1 electron 
system, Fi is the resonance width of the 1 electron system, 
and T2 is the width of the 2 electron system) 



the finite lifetime of the (J ma x + 1)— metastable state, 
the magnitude of this jump is larger than what could be 
expected by setting A — 0. 
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